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Abstract 

Photon trapping and outflow are two key physics associated with the supercritical 
accretion flow. We investigate the conflict between these two processes based on two- 
dimensional radiation-hydrodynamic (RHD) simulation data and construct a simplified 
(radially) one-dimensional model. Mass loss due to outflow, which is not considered in 
the slim-disk model, will reduce surface density of the flow, and if very significant, it will 
totally suppress photon trapping effects. If the photon trapping is very significant, con- 
versely, outflow will be suppressed because radiation pressure force will be reduced. To 
see what actually occurs, we examine the RHD simulation data and evaluate the accretion 
rate and outflow rate as functions of radius. We find that the former monotonically de- 
creases, while the latter increases, as the radius decreases. However, the former is kept 
constant at small radii, inside several Schwarzschild radii, since the outflow is suppressed 
by the photon trapping effects. To understand the conflict between the photon trapping 
and outflow in a simpler way, we model the radial distribution of the accretion rate from 
the simulation data and build up a new (radially) one-dimensional model, which is similar 
to the slim-disk model but incorporates the mass loss effects due to the outflow. We find 
that the surface density (and, hence, the optical depth) is much reduced even inside the 
trapping radius, compared with the case without outflow, whereas the effective temperature 
distribution hardly changes. That is, the emergent spectra do not sensitively depend on the 
amount of mass outflow. We conclude that the slim-disk approach is valid for interpreting 
observations, even if the outflow is taken into account. The observational implications of 
our findings are briefly discussed in relation to ultra-luminous X-ray sources. 
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radiation mechanisms: general — stars: winds, outflows 
1. Introduction 

Supercritical (or super-Eddington) accretion onto black holes remains one of the most funda- 
mental, classical issues in the present-day astrophysics and is now discussed in wide fields of astro- 
physics (see Chapter 10 of Kato et al. 2008 for a concise review). It is well known for the case of 
spherical accretion that there is an upper limit to the luminosity; that is the Eddington luminosity, 
Le- It is thus impossible for gas to accrete onto a black hole at a rate exceeding the critical accre- 
tion rate. Merit [= -^e/ ('7c^)] where c is the speed of light, r] is the efficiency. Then, how about the 
cases of disk accretion? Is the supercritical accretion (accretion at a rate exceeding the Eddington 
rate) feasible? This is the enigmatic issue and has been discussed from the 1970's by many authors, 
including Shakura & Sunyaev (1973), but it still remains as a controversial issue because of technical 
difficulties by the analytical approach. One of the reasons for these technical difficulties stems from 
the multi-dimensional properties of the supercritical accretion flow. 

There are two key processes which appear when the disk luminosity approaches the Eddington 
luminosity: photon trapping and radiation pressure-driven outflow, and both are multi-dimensional 
effects. At very high luminosity, the accretion rate should also be very high, and so is the optical 
depth. Then, the photon diffusion timescale in the vertical direction may become shorter than the 
accretion time of gas. If this happens, photons generated deep inside the accretion flow are not able to 
reach the surface before the material is swallowed by a black hole. This is the photon trapping effects 
(Katz 1977; Begelman 1978; Begelman & Meier 1982; Flammang 1984; Blondin 1986; Colpi 1988; 
Wang & Zhou 1999; Ohsuga et al. 2002). Furthermore, the flow of high luminosities are supported 
by radiation pressure, which is likely to induce outflow (Bisnovatyi-Kogan & Blinnikov 1977; Meier 
1979; Icke 1980; Tajima & Fukue 1998). For complete understanding of the supercritical accretion 
flow, we need to solve the multi-dimensional radiation-hydrodynamic (RHD) equations. This has 
become possible quite recently thanks to the rapid developments of high-speed computers. 

Ohsuga et al. (2005) were the first to succeed in the global RHD simulations of the supercritical 
accretion flow until the flow settles down on the quasi-steady phase. They have demonstrated that the 
accretion rate can be arbitrarily high, that the Eddington luminosity can be exceeded, and that the 
luminosity increases with an increase of the mass input rate in a logarithmic fashion (see also Ohsuga 
2006). The reason for the occurrence of super-Eddington luminosity is the combination of two effects: 
significant radiation anisotropy and photon trapping (Ohsuga & Mineshige 2007). 

The slim-disk model was proposed by Abramowicz et al. (1988) as a simplified model of 
supercritical accretion flow. This model is constructed on (radially) the one-dimensional formulation, 
like Shakura & Sunyaev (1973), and the photon trapping effects are expressed as the (radial) advection 
of radiation entropy, though the outflow effect is not considered. Despite this weakness this model has 
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been used by many authors as a "standard model" of the supercritical accretion because of its technical 
simplicity (e.g., Szuszkiewicz et al. 1996; Beloborodov 1998; Watarai & Fukue 1999; Mineshige et 
al. 2000; Fukue 2000; Watarai et al. 2000; Kawaguchi 2003; Gu & Lu 2007). 

It is not easy to construct one-dimensional models (like the slim-disk model) which incorpo- 
rate both of photon trapping and outflow effects, since the combined effects are essentially multi- 
dimensional. There have been several attempts (e.g., Lipunova 1999; Kitabatake et al. 2002; Fukue 
2004; Kohri et al. 2007; Poutanen et al. 2007), but these studies rely on simplifications and assump- 
tions. It will be much more preferable to study this issue by using the multi-dimensional RHD simu- 
lations, but such simulations are expensive and it is not easy to gain physical insight from simulation 
data. In the present study we thus construct a one-dimensional model as an extension of the slim-disk 
model based on the information regarding the conflict between photon trapping and outflow, which 
we gain from the simulation data. Our goal is to examine how the flow structure and its appearance 
are affected by the counteracting effects of photon trapping and outflow. 

The plan of this paper is as follows: In the next section, we overview the RHD simulations of 
supercritical accretion flow and calculate the accretion rate and the mass outflow rate of the simulated 
flow. In section 3, we present our improved slim-disk model and see some details of the flow structure. 
We then discuss several important effects of the supercritical flow based on our simple model. The 
final section is devoted to summary. 

2. Supercritical Accretion Flows 

2.1. Basic Considerations 

Through the studies based on the slim-disk model two key observational signatures of the 
supercritical flow have become clear: smaller innermost radius and the flatter temperature profile 
(Watarai et al. 2000). The former is because of large amount of accreting material existing even 
inside the radius of the innermost stable circular orbit (ISCO), emitting significant radiation (see also 
Watarai & Mineshige 2003). The latter is due to the suppression of the radiation flux by photon 
trapping. 

The emergent flux distribution of the standard-type disks is determined by the energy balance, 
which is approximately expressed by 

o d2 rpi GMMs.cc n_i ... 

211 R ■ crTgg oc oc R . (1) 

R 

This gives the relation, T^s oc R~^^^. Here, a, G, M, and M^^c represent, respectively, the Stefan- 
Boltzmann constant, the gravitational constant, mass of the central black hole, and the mass accretion 
rate. 

In the slim-disk model, by contrast, the advection effects, the relative importance of which 
increase inward as tdifr/^acc oc R^^ (see Kato et al. 2008, Chap. 10), should be considered. As a 
result, the flux distribution becomes 
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o d2 rr4 GMMg^cc ^acc j-,0 /o\ 

211 R ■ crT^Q OC oc R , (2) 

R ^diflf 

which lead to a somewhat flatter temperature profile, Teg oc Here, tacc (= —R/vr) and tjifr (= 

2>Ht /c) represent the accretion timescale and photon diffusion timescale (in the vertical direction), 
respectively, and vr, H, and r are radial velocity, scale-height of the disk, and optical thickness of the 
disk, respectively. The occurrence of photon trapping gives rise to a critical radius, trapping radius 
(-Rtrap)^ inside which photon trapping is significant, and it is expressed as 

i?trap~^(MacccVLE)rs, (3) 
ix 

where (= 2GM/c^) is the Schewarzschild radius. 

Cautions should be taken to the fact that the slim-disk model does not perfectly describe the 
properties of the supercritical accretion flow, since outflow, one of the most important properties of 
the supercritical flow, is not considered. The mass accretion rate is, hence, set to be constant in 
space in the slim-disk formulation. However, it is likely to decrease inward due to the mass loss 
by the outflow. In the limit of significant mass loss by the outflow, surface density (and thus the 
optical depth) of the accretion flow will be much reduced, which may result in the total suppression 
of the photon trapping effects. That is, the slim-disk model breaks down. In the limit of significant 
photon trapping, conversely, (outward) radiation pressure force will be weakened (or may become 
even inward, see Ohsuga & Mineshige 2007), and the mass loss by the radiation pressure-driven 
outflow will be negligible. Which is likely to be the case? Which process will dominate? 

We wish to note that the slim-disk model is not the only models of the supercritical accretion. 
Shakura & Sunyaev (1973) already discussed the supercritical flow and considered standard-type 
disks with large outflow. They assumed that flux from each radius F{R) cannot exceed the value 
that gives the Eddington luminosity, i.e., roughly 2'kR^F{R) = within certain radius, called the 
spherization radius, Rgpt (see also Begelman 1979). Using the standard-disk relation [Eq. (1)], we 
can easily derive 

i?sph~ (MacccV^E)rs (4) 

Inside this radius, the mass accretion rate decreases inward in proportion to the radius; i.e.. Mace oc R. 
This means, the accretion rate vanishes at a very small R, and we again have a somewhat flatter 
temperature profile, Tcs oc since 27ri?^ ■ aT^g oc Ms.cc/ R^ const. 

Comparing equations (3) and (4), one can notice that the critical radii for the photon trapping 
and outflow are on the same order because H R. In other words, both effects could be equally 
important in the supercritical flow. In addition, the apparent effective temperature profiles are same 
for both cases. That is, we cannot simply conclude which effect dominates by looking at the observed 
spectra. 

Theoretically speaking, however, there is a big difference. The photon trapping is more im- 
portant near the equatorial plane, since it takes longer time for photons to travel from the equatorial 
plane, than from the middle part of the disk, to the disk surface, while outflow occurs from the disk 
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surface. This indicates that both effects can be simultaneously significant but at different heights. We 
thus need to carefully examine multi-dimensional simulation data. 

2.2. Model and Simulated Flow 

We simulate the supercritical accretion flow, by using the two-dimensional RHD code devel- 
oped by Ohsuga et al. (2005). The basic equations and the numerical method are described in details 
in Ohsuga et al. (2005). Hence we briefly summarize them. We use spherical coordinates (r,9,^p), 
where r, 9, and if are the radial distance, the polar angle, and the azimuthal angle, respectively, and 
set a non-rotating black hole at the origin. To mimic the general relativistic effects, we adopt the 
pseudo-Newtonian potential, ^ given hy ^ = —GM/ (r — Ts) (Paczyhski & Wiita 1980). As to the 
flow structure we assume the axisymmetry (i.e., d/d(^ = 0) and the reflection symmetry relative to 
the equatorial plan (with 6 = n/2). To solve the radiative transfer, we apply a flux-limited diffusion 
(FLD) approximation developed by Levermore & Pomraning (1981). This approximation is that the 
radiative flux and the radiation pressure tensor are expressed in terms of the radiation energy density 
(Turner & Stone 2001; Ohsuga et al. 2005). 

A difference between the present calculation and Ohsuga et al. (2005) is in the computational 
domain. Our purpose in this study is to examine the conflict between outflow and photon trapping. 
Therefore, we have to simulate the supercritical accretion flow in a wider spatial range. Thus, we set 
the computational domain spherical shells of 3rs < r < 1000 rs and < 9 < n/2 and it is divided into 
96 X 96 grid cells, (Note that the computational domain of Ohsuga et al. 2005 was 3rs < r < 500 rg.) 
We start the calculations with a hot, rarefied, and optically thin atmosphere. There is no cold dense 
disk initially, and we assume steady mass injection into the computational domain through the outer 
disk boundary (r = lOOOrg, 0.457r < 9 < O.Stt). The injected matter is supposed to have a specific 
angular momentum corresponding to the Keplerian angular momentum at r = 500 rs (cf. r = 100 rs in 
Ohsuga et al. 2005), and we set the injected mass-accretion rate (mass input rate) Minput to be constant 
in time. Throughout the present simulation, we assume M = IOM0, a = 0.1,7 = 5/3,/i = 0.5, and 
Z = Zq. Here a, 7, /i, and Z are the viscous parameter, the specific heat ratio, the mean molecular 
weight, and the metallicity, respectively. 

Figure 1 indicates the color contours of the matter density (left panel) and the radial inflow 
velocity (right panel) distribution in the meridional plane, which are time-averaged over t = 190 — 
250 [s] in the case of Mmput = 1000 Le/c^. In the right panel the radial inflow velocity is normalized 
by the escape velocity and the region with white color indicates the outflow region, i.e., Vr > 0. The 
supercritical accretion flow forms at r < 350 rg. The disk accretion of the high density gas and the 
outflow of the low density one are clear. In the previous calculation, the supercritical accretion flow 
forms at r ^ SOrg. The behavior of the present simulated flow in this region is roughly consistent with 
that of the previous one (Ohsuga et al. 2005; see also Ohsuga & Mineshige 2007). 
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Fig. 1. Color contours of the matter density {left) and the radial inflow velocity {right) distribution in the meridional 
section. Note that these values are time-averaged over t = 190 — 250 [s]. 




Fig. 2. Schematic picture for the calculation of mass-accretion/outflow rate. 

2.3. Mass Accretion/Outflow Rate 

We assume that an accretion flow occurs in the region of 6'disk < 6* < 7r/2, while outflow region 
corresponds to the region above the accretion flow region, and determine an angle 6'disk (> 0) from the 
simulation data. That is, 6'disk is chosen so as to roughly coincide with the angle between the 2;-axis 
and the boundary separating the region with radial inflow (negative velocity) and that of the positive 
velocity (figure 2). On the basis of this assumption, we calculated the mass accretion rate and the 
cumulative mass outflow rate as functions of radius. Since the quasi-steady flows form at r < 350 rg, 
we set in this calculation the outer radius of the flows to be r„ind = 350 Ts. 

We then calculate the accretion rate, Macc('^), and the cumulative mass outflow rate, Mout('^)> 

by 

/■90° 

Macc(r)=/ 47rrVrSin6'rf6', (5) 

and 

Mont{r)= I 47rr pt^e 1 0^.^^ sin 6'disk . (6) 

J r 
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In figure 3, we show the accretion rate and mass outflow rate as functions of radius in the sim- 
ulated flows for Minput = 1000Le/c^ {left panel) and Mmput = 3000Le/c^ {right panel), respectively. 
These values are normalized by the Eddington accretion rate, i.e., L-e/c^. Here, the disk inclination 
angle is chosen to be ^disk = 70° in the case of Minput = 1000 Le/c^, and 6'disk = 60° in the case of 
Minput = 3000 Le/c^. In each panel, the solid, dashed, and dotted curves represent the accretion rate, 
Macc('"), the cumulative mass outflow rate, Mont{f) , and the sum of the these rates, respectively. In 
particular, the following mass conservation should hold. 

Mace (rwind) = Mace (r) + Mout (r) = const. (7) 

We confirm this relation within errors less than 4%. The reason why Maec(r„ind) < Minput is because 
a part of injected gases accumulates at r^ind < r < 1000 rg or goes out to the outer region of the 
computational domain without accreting immediately. 

The accretion rate is not constant in space due to the mass loss by outflows. In the analytical 
accretion flow model, in which an outflow is considered, the accretion rate decreases inward in pro- 
portion to the radius, i.e., M^ce oc r (e.g., Shakura & Sunyaev 1973). That is, the accretion rate should 
vanish at a very small radius. However, we find that the accretion rates still have a finite value at small 
radii, inside several Schwarzschild radii. This is because the emergence of outflow is suppressed due 
to the attenuation of radiation flux by photon trapping. In fact, the radiation flux become even neg- 
ative (inward) at these radii (Ohsuga & Mineshige 2007). The outflows are difficult to blow off and 
thus the flow is easy to be accreted towards the central black hole. 

In the case of A/input = 1000 Le/c^, the outflow blows off at r ^ SOOrg. As we evaluated the 
photon trapping radius by equation (3), we found that the photon trapping is effective at r < 150 rg 
for a mass accretion rate of Mace(r„ind = 350rs) = 364Le/c^. In the case of Minput = 3000Le/c^, in 
contrast, the outflow and photon trapping is already effective at r„ind = 350rs. Thus, we confirmed that 
the critical radii for the photon trapping and outflow are on the same order in the RHD simulations. 
Note that this is only a rough estimate, since we assumed constant accretion rate in the derivation. 

3. Modified Slim-Disk Model 

On the basis of the simulation data analysis presented in the previous section, we try to con- 
struct a modified slim-disk model, which incorporates the effects of the mass loss by an outflow, in 
this section. The accretion flow models which incorporate the mass loss effect were proposed by 
several researchers (Lipunova 1999; Kitabatake et al. 2002; Fukue 2004; Poutanen et al. 2007). They 
evaluate the outflow by adopting the spherization radius. By contrast, we construct a model by using 
the accretion rate obtained in the previous section. 

3.1. Basic Equations 

We use cylindrical coordinates (-R, v^, z) in this section. We assume steady and axisymmetric 
flow, and non-rotating black hole, and adopt a pseudo-Newtonian potential. We use height-integrated 
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Fig. 3. The accretion rate and cumulative mass outflow rate as functions of radius in the simulated flows in the case 
of A/input = IOOOXe/c^ {leji) and Minput = 3000Le/c-^ {right) . The sohd, dashed, and dotted curves represent the 
accretion rate Afacc(?'), mass outflow rate Mout(?') , and the sum of the these rates, respectively. The data of both 
simulated flows is time-averaged over t = 190 — 250 [s]. The accretion rates at the outer radius rwind (= 350rs) are 
A^acc(?'wind) = 364iyE/c^ {l^ft) and 828iE/c^ {right). The shadowed areas indicate the photon trapping region. 

quantities, such as S = J^jjpdz = 2I^pH and 11 = J^jjpdz = 2I^j^ipH. Here, SAIlAp, and H 
are the surface density , height-integrated pressure, total pressure, and scale height, respectively. The 
coefficient Jn and Jn+i were introduced by Hoshi (1977). The density and the pressure are related to 
each other by the polytropic relation, p oc p^+^Z^. We assign N = 3 throughout the entire calculation 
(i.e., Js = 16/35 and h = 128/315). 

The continuity equation, the radial component of the momentum equation, the angular mo- 
mentum conservation, the hydrostatic equilibrium in the vertical direction, the energy equation, and 
equation of state are written as follows: 

Mace (R) = Mace (R^ind) ' Mout (R) , (8) 

dvR IdU n rflnl^K 

^^dR ^ S Zr ~ ~R3 S dR ' ^ ^ 

M^cc{i-iin) = -27rR'Tn^, (10) 

{2N + 3)^ = H'nl, (11) 

Qtis = Qvad + Qadv (12) 

n = Hgas + Hrad = ^^ST^ + \h^^aT'^H. (13) 

Here, the flow velocity of radial component and azimuthal one are expressed by vr and v^, respec- 
tively, the angular momentum of the gas is given by £(= Rv^ = R^^l), and I^k are the angular 
speed of rotation and the Keplerian angular speed, Ek and £in are Keplar angular momentum and the 
angular momentum at inner radius of the flow, Tr^{= —all) is vertically integrated stress tensor with 
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a being the viscosity parameter (Shakura & Sunyaev 1973), Ilgas and Ilrad are the gas pressure and 
the radiation pressure, and niu, Tc, ks, and a are the hydrogen mass, the temperature on the equa- 
torial plane, the Boltzmann constant, and the radiation constant, respectively. The last term on the 
right-hand side of equation (9) is a correction term resulting from the fact that the radial component 
of the gravitational force changes with height (Matsumoto et al. 1984). 

In the energy equation (12), the viscous heating rate Q^^^, the radiative cooling rate Q~d' ^^'^ 
the advective cooling rate Qadv defined by 

Q+^ = -RTn^^, (14) 
a-ad = 2F=^, (15) 



g^dv = o^'^s^c— , (16) 



9 ds 
'dR' 

where F is the radiative flux per unit surface area on the flow surface, s is the specific entropy, and 
r is the optical thickness of the flow given by r = {k^s + f^s)^- 'tcs(= 0.4) is the electron scattering 
opacity, Kfj(= 0.64 x 10^^ pT"''/^) is the free-free absorption opacity, and p and T are the vertically 
averaged density and the temperature, respectively. Note that Jiao et al. (2009) point out that the usage 
of equation (11) may lead to overestimation of gravity force when the scaleheight is comparable to 
the radius, H ^ R. 

We wish to note that the difference between the equations of the original slim-disk model 
and those of our model is only in the continuity equation (8). In the original slim-disk model, the 
continuity equation is expressed by Macc(-R) = const. The reason why the forms of equations (9)- 
(13) do not change is because these expressions are written per unit mass. While Lipunova (1999) 
explicitly expresses the outflow effect in her equations, the physical meanings of the equations are 
the same as own (see also Kitabatake et al. 2002). The difference between the previous models and 
our model is that we construct a model which realistically considers the outflow effect by using the 
accretion rate in the previous section. In the previous study, moreover, the flow is approximated to be 
Keplarian (Lipunova 1999; see also Poutanen et al. 2007) or assumed to obey a self-similar solution 
(Kitabatake et al. 2002; see also Fukue 2004), whereas we solve numerically the radial advection, i.e., 
equation (9) and (12). However, we need to take into account the work exerted on the outflow by the 
accretion flow material. This issue will be discussed in the next section. 

The calculations were performed from the outer radius at i? = 1 x 10^ rg down to the in- 
ner radius, R ~ 1.0 rg. In our modified slim-disk model, we approximate M^cdR) = M^cdr) and 
Mout{R) = Mout{r) within R = 350 rg. Since the range of the accretion rate which is obtained 
by the simulation in the previous section is at i? < 350 rg, we assume that the accretion rate at 
R = 350 rs - 10 is constant, i.e., M,,,{R') = Macc(i?wind) = M,cc{r^ind) for 350 < R' < 10^ r,. 
In the original slim-disk model, in contrast, Macc(-R) = Afacc(-Rwind) in the whole radius. We set the 
black hole mass and the viscous parameter to be M = IOM0 and a = 0.1, respectively. We also set 
Macc(i?wind) = 364Le/c2 and 828 Le/c\ 

9 





Fig. 4. The surface density E and the effective temperature Tcs as functions of radius in the modified 
slim-disk model (dashed lines) and the original slim-disk model (solid lines) for A/input = 1000 Le/c^ (left) and 
A^i„put = 30G0iE/c2 {right). 

3.2. Flow Structure 

Figure 4 indicate the surface density and the effective temperature as functions of radius in the 
modified slim-disk model (dashed lines) and the original slim-disk model (solid lines) for Minput = 
1000 Le/c^ (left panels) and Minput = 3000 Le/c^ (right panels), respectively. The surface density of 
the flow is significantly reduced due to the mass loss by outflow in both cases (see figure 3). Hence, 
the optical depth of the flow r (oc S) also is reduced by a factor of two or three. However, the flow 
is optically thick in the whole region. Even though the mass accretion rate decreases as the radius 
decreases, the trapping radius derived from the mass accretion rate at any radius exceeds that radius; 
i.e., photon trapping is effective. 

We quantitatively examine the effect of photon trapping, i.e., we calculate the relative impor- 
tance of the advective cooling in the energy equation. We derive the relation between tdifr/^acc and 
Qadv/Qvis- The advective cooling rate are expressed by 
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Fig. 5. Ratio of the advective cooling rate Q^dv '■° '■^^ viscous heating rate Q^^^ as functions of radius in the 
modified sUm-disk model (dashed lines) and the original sUm-disk model (solid lines) for A/jnput = 1000 Le/c^ 
(left) and A/input = 3000 Le/c^ (right) . 

Qadv--8^S(^e--j=---i/aT,, (17) 

where e is the specific energy. Here, we used e = aT^/p and p = aT^/3. Therefore, the condition of 
effective photon trapping (tdiir/^acc > 1) is derived by 

tdiff 32 Qadv 32 Qadv ^ t ^1 q\ 

''acc -J Vrad ^ Vvis Vadv 

As a result, this gives the relation, 

^> — -0.09. (19) 

Qvis 35 

That is, photon trapping is effective at Q^^/Qtis ^ 0.09. 

In figure 5, we show the ratio, Q^dv/Qtis^ each model. In the case of the original slim-disk 
model, in which the outflow is not considered, photon trapping is significantly effective at i? ^ 100 rg. 
In the case of the modified slim-disk model, however, the advective cooling rate is smaller due to the 
mass loss by outflow (Qadv '^he ratio. Hence, the photon trapping effect is weaker in 

the model, but the outflow is not strong enough to totally suppress photon trapping at the smaller radii 
(Qadv/Qtis ~ 0-4 — 0.6). At the larger radius, there is no difference between the two models because 
there is no outflow. 

In addition, we understand from figure 4 that, even if outflow effects are taken into account, the 
effective temperature profiles do not change. It is Tcff oc r~^/^ at the smaller radius, and Tcs oc r~^/^ 
at the larger radius. Why is then the effective temperature profile unchanged? This is because it is 



determined by the photon trapping at small radii (note that photon trapping is effective at Qadv/ Q'vis ~ 



0.09), while no significant outflow occurs at large radii . 

Kitabatake et al. (2002) construct a model for supercritical accretion flows with mass loss, 
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adopting the self-similar treatment proposed by Narayan & Yi (1994). They mention that the effective 
temperature profile of the flow is negligibly affected by occurrence of outflow in the case of advective- 
dominated flows (see also Fukue 2004). The reason is as follows. When the radiative flux F (= aTc oc 
IT) is reduced by the mass loss, optical depth r (oc S) is reduced at the same time. As a result, the 
effective temperature of the flow Tcs (= TcT"^/"') does not depend on the mass loss. This provides 
another explanation for the flatter temperature distribution. 

In the region at i? ^ 3rs, the advective cooling rate becomes negative. The reason is as follows. 
The optical depth of the flow r (oc S) steeply decreases inward in that region, although the flow is 
optically thick (see figure 4). Therefore, the radiative cooling rate (oc r"^) steeply increases 
inward, whereas the viscous heating rate (Qtis ^^acc) does not. As a result, the advective cooling 
rate should become negative, since Qadv = Qtis ~ Q7axi- 

4. Discussion 

4.1. Brief summary 

In this paper, we first carefully examined the global RHD simulation data of supercritical 
accretion flow onto black holes in order to examine the conflict between photon trapping and outflow. 
We have confirmed that both are equally important; i.e., despite significant mass loss by the outflow, 
it is not strong enough to totally suppress photon trapping. We evaluated the accretion and outflow 
rates as functions of radius based on the simulation data and put them into the formulation of the slim- 
disk model, in which no outflows were considered originally and hence the mass accretion rate was 
considered to be constant. We compared the resultant flow structure with consideration of outflow 
by those without consideration of outflow, finding that, although surface density (and, hence, optical 
depth) of the flow is significantly reduced, the effective temperature profile is negligibly affected 
by the occurrence of outflow. Therefore, multi-blackbody spectra are negligibly affected, either. 
This has a profound implication when one performs spectral fitting of black hole objects, notably 
of Ultraluminous X-ray Sources (ULXs, see subsection 4.3). 

We here remark on the reasons why we stick to the one-dimensional model, when multi- 
dimensional simulation data are available. There are a number of reasons for this. Although it has 
become possible to simulate the flow from the first principle, they are still subject to numerical errors 
and limitations arising from the finite mesh spacing. Also, it is not always easy to specify the physical 
processes from the vast simulation data with substantial fluctuations and numerical errors. Further, 
multi-dimensional RHD simulations are very expensive and time-consuming. That is, it is impossible 
to perform extensive parameter studies. Therefore, studies based on simplified (one-dimensional) 
models like the present one should be useful and beneficial for understanding the physics. 

4.2. Work exerted on outflow 

In our simplified one-dimensional model, we consider loss of mass, angular momentum, and 
energy by outflow. These effects are incorporated by spatial variation of the mass accretion rate. I.e., 
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we considered angular momentum and energy loss carried by outflow material, assuming specific 
angular momentum and specific energy are the same for the outflow and the accretion flow. There 
exists, however, another important factor which should be include in the energy equation; that is the 
work exerted on the outflow by the disk material, Q~ind- Obviously, outflow cannot have positive 
energy to reach the infinity without acquiring additional energy by the underlying accretion flow. 
Poutanen et al. (2007), for example, include this effect by assuming that radiative loss of the accreting 
material is partly transferred to the outflow. 

To see how our results are affected by this additional energy loss, we reduce the specific energy 
of accretion flow by hand, assuming Q^i^^ = Q^ad- The results are that the effective temperature 
profile hardly varies, even though the amount of mass loss decreases. Our conclusion of unchanged 
temperature profile is not altered. 

4.3. Model for Ultraluminous X-Ray Sources 

ULXs are bright, compact X-ray sources found in the off-center region of nearby galaxies. 
Their luminosities are typically, Lx ~ lo^^"^^ [erg/ s], and, hence, exceed the Eddington luminosity of 
a neutron star (Fabbiano 1989). There are two hypotheses for the origin of ULXs. If the luminosity 
is below the Eddington luminosity, it then follows that the mass of the central black holes of ULXs 
should be intermediate-mass black hole (IMBH), whose masses range over lO^'^M©. If the lumi- 
nosity can be above the Eddington luminosity, conversely, stellar-mass black holes can also account 
for the luminosities. Unfortunately, the kinematic method, which is very powerful to estimate the 
black hole masses in binary systems, is not applied to the ULXs because no appreciable (optical) line 
features are observed. 

Interestingly, ULXs share similar spectral features with the black hole binaries (Colbert & 
Mushotzky 1999; Makishima et al. 2000; see review by Done et al. 2007). Therefore, it seems possible 
to estimate the black hole mass through the spectral fitting, although the results look controversial. 
Some authors claim that the ULXs should have IMBHs. This is because the obtained blackbody 
temperature is low, on the order of ~ 0.1 keV (Cropper et al. 2004; Miller et al. 2004; Roberts et al. 
2005). (The blackbody temperature is proportional to M~^/^ for the same Eddington ratio.) Some 
others claimed that the ULXs should contain stellar-mass black holes (King et al. 2001; Watarai et 
al. 2001; Okajima et al. 2006; Vierdayanti et al. 2006). Vierdayanti et al. (2008), for example, claim 
that the spectral fitting with the conventional spectral model (with blackbody and power law spectral 
components) is not reliable when the power-law component dominates, and demonstrated basing on 
the original slim-disk model that some ULXs exhibit spectral signatures of the supercritical accretion 
flow (see also Vierdayanti et al. 2006). The present study supports their conclusions, since, even if 
outflow effects are taken into account, neither of the effective temperature profiles nor multi-color 
blackbody spectra are altered. 
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